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ABSTRACT 



The kinematics of 129 W UMa binaries is studied and its implications on the 
contact binary evolution is discussed. The sample is found to be heterogeneous in 
the velocity space that kinematically younger and older contact binaries exist in the 
sample. Kinematically young (0.5 Gyr) sub-sample (MG) is formed by selecting the 
systems which are satisfying the kinematical criteria of moving groups. After removing 
the possible MG members and the systems which are known to be members of open 
' clusters, the rest of the sample is called Field Contact Binaries (FCB). The FCB 

, has further divided into four groups according to the orbital period ranges. Then a 

correlation has been found in the sense that shorter period less massive systems have 
larger velocity dispersions than the longer period more massive systems. Dispersions in 
the velocity space indicates 5.47 Gyr kinematical age for the FCB group. Comparing 
■ with the field chromospherically active binaries (CAB), presumably detached binary 

progenitors of the contact systems, the FCB appears to be 1.61 Gyr older. Assuming 
O ■ an equilibrium in the formation and destruction of CAB and W UMa systems in the 

I 1 Galaxy, this age difference is treated as empirically deduced lifetime of the contact 

stage. Since the kinematical ages (3.21, 3.51, 7.14 and 8.89 Gyr) of the four sub 
groups of FCB are much longer than the 1.61 Gyr lifetime of the contact stage, the 
pre-contact stages of FCB must dominantly be producing the large dispersions. The 
kinematically young (0.5 Gyr) MG group covers the same total mass, period and 
spectral ranges as the FCB. But, the very young age of this group does not leave enough 
room for pre-contact stages, thus it is most likely that those systems were formed in the 
beginning of the main-sequence or during the pre-main-sequence contraction phase, 
either by a fission process (Roxburgh 1966) or most probably by fast spiraling in of 
two components in a common envelope. 

Key words: binaries: eclipsing, stars: binaries spectroscopic, stars: evolution, stars: 
kinematics 



1 INTRODUCTION 

The low mass contact binaries, popularly known as W Ursae 
Majoris stars, are easily recognized by their eclipsing light 
curves with equal (or nearly equal) minima as wide as to 
touch one another. With a small orbit implied by a short 
orbital period less than a day, the binary components fill 
their Roche lobes so that the dumbbell-shaped stars touch 
each other at the inner Langrangian point. The W UMa 
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systems have spectral types ranging from late A to mid- 
dle K and they are located near or just above the main 
sequence. The strong tidal interaction synchronized their ro- 
tation thus the period of rotation is equal to the orbital rev- 
olution. The rapid rotation and convective atmospheres are 
primary causes to their observed chromospheric and coro- 
nal emissions as well as to the presence of large starspots 
which are the manifestations of strong, dynamo-generated 
magnetic activity. 

The rapid rotation and quickly changing radial velocity, 
on the other hand, is a serious obstacle to obtain a reliable 
radial velocity from the highly broadened and blended spec- 
tral lines of numerous W UMa systems. Therefore, forming 
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a complete sample for studying galactic space motions and 
kinematics of W UMa stars is limited with the availability 
of the radial velocity measurements. From the identified 374 
W UMa binaries (Kukarkin et al. 1971; van't Veer 1975) in 
the 3rd General Catalog of Variable Stars (GCVS), Guinan 
& Bradstreet (1988) were able to collect only 34 systems 
with a measured center of mass velocity. The number of 
stars classified as EW (eclipsing binaries of W UMa type) 
increased to 561 in the 4th edition of GCVS. But there were 
only 78 systems available with complete physical parame- 
ters which were studied statistically by Maceroni & van't 
Veer (1996). Only 37 out 78 had spectroscopic mass ra- 
tio determined from the radial velocities. Finally, Asian et 
al. (1999) counted 125 W UMa binaries in the Hipparcos 
Catalogue, but there were only 34 binaries from which the 
space motions can be computed. The revised 4th edition 
of GCVS now contains 751 W UMa. Thanks to the dedi- 
cation of David Dunlap Observatory (DDO), where 54 per 
cent of current radial velocities came from, and Dominion 
Astrophysical Observatory (DAO) with 20 per cent contri- 
bution, and the improvements of spectroscopic techniques 
for measuring radial velocities that a total of 129 W UMa 
systems were collected for this study. Nearly four times the 
number and the improved accuracy of data justified us to 
repeat previous attempts of studying kinematics of W UMa 
binaries. 

With a small sample of less accurate data, Guinan & 
Bradstreet (1988) estimated a kinematical age of 8-10 Gyr 
for W UMa systems. If this age is compared to the ~ 5 
Gyr average age of chromospherically active binaries (CAB), 
classical RS CVn and BY Dra systems (Eker 1992), W UMa 
stars appear two times older. Nevertheless, this age of W 
UMa binaries was found consistent with the theory of con- 
tact binary formation from the detached binary progenitors, 
presumably from the short period (P < 5 days) RS CVn bi- 
naries, by a mechanism involving angular momentum loss 
and orbital period decrease previously suggested by Huang 
(1967), Mestel (1968), Okamoto & Sato (1970), van't Veer 
(1979) and Vilhu & Raunen (1980). But, according to the 
space velocities and dispersions computed by Asian et al. 
(1999), W UMa binaries cannot be older than RS CVn 
systems. This means that W UMa binaries could not have 
evolved from detached binaries unless the time to reach con- 
tact stage is very short. 

With no pre-contact stage, there is an opposing theory 
which predicts the formation of contact binaries directly by 
a fission process (Roxburgh 1966) at the end of the pre- 
main sequence contraction. Considering the brief lifetime of 
the contact stage (0.1 < t con tact < 1 Gyr) estimated by 
Guinan & Bradstreet (1988), this theory would fail to pro- 
duce older populations of W UMa binaries unless the an- 
gular momentum of contact binaries is conserved. However, 
with the conservation of the orbital angular momentum, this 
theory would refuse the formation of the contact binaries 
from the detached binary progenitors. 

Apparently, an empirical evidence implying whether the 
total orbital angular momentum is conserved or not within 
the evolutionary time scale is extremely important to re- 
veal the actual working model. Inspired by the predictions 
of Demircan (1999), Karatas, et al. (2004) presented such an 
evidence recently from the kinematics of CAB sample con- 
taining RS CVn and BY Dra stars. The evidence of orbital 



period decrease, connected to mass and angular momentum 
loss, was shown by comparing the total mass (Mi + M2) 
and period histograms between the kinematically younger 
(age 0.95 Gyr) and the older (age 3.86 Gyr) CAB sub sam- 
ples. The orbital period decrease was also indicated by the 
fact that the shorter orbital period systems were found older 
than the longer orbital period systems. According to a pre- 
liminary investigation of Demircan et al. (2004), the period 
decrease is exponential in time so that the active detached 
binaries, up to about 5 days orbital period, could change 
into contact form within less than about 5 Gyr. Therefore, 
the formation of W UMa stars from detached binaries is ev- 
idently occurring. Nevertheless, this does not mean this is 
the only mechanism to form contact binaries. It may still 
possible that a limited fraction W UMa systems, possibly 
the very young ones, may be forming by the fission process. 

In the present work, the kinematics of W UMa bina- 
ries is studied and its implications on the contact binary 
evolution is discussed. The present sample has been found 
heterogeneous like the CAB sample in a sense that there are 
younger (~ 1 Gyr age) and older (~ 6 Gyr age) W UMa bi- 
naries covering the same range of spectral types and orbital 
periods. The evidence of mass and angular momentum loss 
causing orbital period decrease during contact binary evo- 
lution has been investigated for a preliminary step as done 
by CAB sample by Karatas, et al. (2004) by comparing the 
total mass (Mi + M2) and period histograms between the 
younger and the older sub samples. 



2 DATA 

Among the 751 recorded W UMa binaries in the recently 
revised 4th edition of GCVS, 129 systems were found to 
have radial velocities. After compiling the rest of the basic 
data (parallaxes and proper motions), the sample is listed 
in Table 1, which is self explanatory, with the columns: or- 
der number, name, HD and Hipparcos reference numbers, 
celestial coordinates (ICRS 2000), proper motions, paral- 
laxes and radial velocities. The associated standard errors 
are given besides the related parameter. The reference num- 
bers in the last column are separated into three fields with 
semicolons to indicate from where the basic data were taken. 
The two or more reference numbers in a field indicate the 
sub fields, that is, there are multiple references. 

2.1 Parallaxes and proper motions 

The parallaxes and the proper motions in Table 1 were 
mostly taken from the Hipparcos and Tycho Catalogs (ESA 
1997) and the Tycho Reference Catalog (Hog et al. 1998). 
The parallax and the proper motion components of three 
stars (TX Cnc, EE Cet, BL Eri) were taken from Roese & 
Bastian (1988). Among the 129 systems in Table 1, the 25 
binaries do not have any trigonometric parallax, and the 12 
binaries have a large relative uncertainty (a^/n > 0.5) in 
the parallax measurements. Considering the fact that Hip- 
parcos measurements are reliable up to 500 pc (two-sigma 
detection limit a ~ 1 mas, Perryman et al. 1997), we de- 
cided to discard the parallaxes of five (AQ Tuc, TY Pub, 
XZ Leo, AG Vir, RZ Com) systems since their measured 
parallaxes indicate a distance further than 500 pc. These 
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Table 1. Hipparcos astromctric and radial velocity data of the W UMa. 
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00 


(4;5;39) 


55 


V752 Ccn 


101799 


57129 


11 


42 


48 


08 


-35 


48 


57.51 


-53 


65 


1 


63 


-24 


28 





93 


9 


51 


1 


47 


29 


10 





82 


(1;1;40) 


56 


AG Vir 


104350 


58605 


12 


01 


03 


50 


+ 13 


00 


30.02 


-2 


46 


1 


46 


-17 


64 





56 


5 


02 





63 


-5 


91 





73 


(1;5;41,58) 


57 


HX UMa 


104425 


58648 


12 


01 


33 


15 


+43 


02 


29.94 


65 


86 


2 


56 


-21 


34 


2 


15 


6 


68 


3 


01 


-19 


88 


1 


11 


(i;i;8) 


58 


CC Com 






12 


12 


06 


20 


+ 22 


31 


58.00 


-120 


00 


2 


50 


19 


30 


2 


50 


12 


49 


1 


56 


-10 


20 


1 


31 


(2;5;42) 


59 


AH Vir 


106400 


59683 


12 


14 


21 


00 


+ 11 


49 


09.40 


47 


74 


3 


03 


-107 


26 


1 


39 


10 


86 


3 


11 


6 


60 





90 


(1;1;43) 


60 


II UMa 


109247 


61327 


12 


32 


54 


84 


+ 54 


47 


42.89 


-42 


50 


1 


57 


-12 


12 


1 


45 


5 


04 


1 


81 


-8 


02 


1 


10 


(i;i;i7) 


61 


RW Com 




61243 


12 


33 


00 


28 


+26 


42 


58.38 


-125 


45 


3 


46 


-36 


05 


1 


89 


11 


45 


2 


45 


-53 


00 


1 


15 


(1;1;44) 


62 


RZ Com 




61414 


12 


35 


05 


06 


+23 


20 


14.02 


17 


44 


1 


79 


-11 


01 


1 


18 


4 


70 





59 


-1 


80 


1 


67 


(1;5;28) 


63 


SX Crv 


110139 


61825 


12 


40 


15 


04 


-18 


48 


00.93 


37 


93 


1 


21 


-5 


63 





92 


10 


90 


1 


21 


8 


71 





94 


(1;1;14) 


64 


BI CVn 




63701 


13 


03 


16 


41 


+36 


37 


00.59 


-6 


38 


1 


23 


-16 


68 


1 


04 


5 


31 


1 


65 


-5 


30 


1 


80 


(1;1;45) 


65 


KZ Vir 


114726 


64433 


13 


12 


22 


87 


+02 


39 


13.78 


-8 


53 





98 





94 





48 


6 


53 





93 


-4 


50 





48 


(1;1;14) 



five stars and ER Ori, which were erroneously listed with 
a negative parallax, are treated like the other 25 systems 
without a trigonometric parallax. For them, the photomet- 
ric parallaxes, which were determined from the period-color- 
brightness relation M(V) = -4.44 log P+3. 02(5-^)0+0. 12 
by Rucinski (2004), were used. The photometric parallaxes 
were calculated with the corrections of the interstellar ab- 
sorption, which were computed from the color excess as 



Ay = 3.1E(B — V). Since all systems were assigned a spec- 
tral type, the color excesses were available together with 
their observed colors. 

The Hipparcos and Tycho catalogues usually give an 
associated uncertainty for all of the parallax and the proper 
motion component measurements. However, there are 18 
systems in the list without an uncertainty at the proper 
motion components. With an optimistic approach, we have 
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4 Bilir et al. 

Table 1 - continued 



ID Name HD HIP a (2000) 4 (2000) ^ a cosS p s 7r 7 References 

(h m s) (° ' ") (mas yr~ 1 ) (mas yr~ 1 ) (mas) (km s~ 1 ) 



66 


KR Com 


115955 


65069 


13 


20 


15 


78 


+ 17 


45 


56 


99 


-6 


70 





78 


-82 


30 





47 


13 


07 





87 


-7 


86 





38 


(1 


1 


17) 


67 


HT Vir 


119931 


67186 


13 


46 


06 


75 


+05 


06 


56 


27 


-107 


79 


3 


24 


-11 


16 


1 


68 


15 


39 


2 


72 


-23 


38 





68 


(1 


1 


35) 


68 


XY Boo 




67431 


13 


49 


11 


57 


+20 


11 


24 


59 


-23 


32 


1 


87 


26 


91 





89 


2 


94 


1 


70 


2 


11 


2 


02 


(1 


1 


28) 


69 


V757 Cen 


120734 


67682 


13 


51 


55 


74 


-36 


37 


24 


82 


-121 


51 


1 


20 


-01 


89 





94 


14 


18 


1 


10 


50 


00 


2 


50 


(1 


1 


46) 


70 


RR Cen 


124689 


69779 


14 


16 


57 


22 


-57 


51 


15 


65 


-52 


69 





77 


-22 


37 





57 


9 


76 





85 
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00 





41 


(1 


1 


30) 


71 


VW Boo 




69826 


14 


17 


26 


03 


+ 12 


34 


03 


45 


5 


45 


2 


09 


-48 


79 





95 


3 


44 


1 


90 


23 


90 





82 


(1 


1 


47) 


72 


NN Vir 


125488 


70020 


14 


19 


37 


74 


+05 


53 


46 


66 


-54 


17 


1 


09 


-19 


63 





60 


9 


48 


1 


14 


-6 


24 





65 


(1 


1 


26) 


73 


EF Boo 


234150 


71107 


14 


32 


30 


54 


+50 


49 


40 


70 


-33 


23 


1 


17 


-37 


34 





93 


6 


00 


1 


06 


-13 


16 


1 


07 


(1 


1 


14) 


74 


CK Boo 


128141 


71319 


14 


35 


03 


76 


+09 


06 


49 


39 


70 


91 


1 


18 


-86 


58 





80 


6 


38 


1 


34 


37 


04 





75 


(1 


1 


26) 


75 


GR Vir 


129903 


72138 


14 


45 


20 


26 


-06 


44 


04 


13 


-84 


94 


1 


47 


47 


34 





99 


18 


83 


1 


18 


-71 


72 





89 


(1 


1 


26) 


76 


AC Boo 




73103 


14 


56 


28 


34 


+46 


21 


44 


14 


-9 


11 


1 


30 


8 


99 


1 


26 


7 


58 


1 


27 


-8 


60 


2 


20 


(1 


1 


36) 


77 


TY Boo 






15 


00 


46 


94 


+35 


07 


54 


89 


-58 


80 


2 


50 


32 


30 


2 


50 


5 


82 





73 


-41 


90 


1 


55 


(2 


5 


48) 


78 


44 Boo 


133640 


73695 


15 


03 


47 


30 


+47 


39 


14 


62 


-436 


26 


1 


36 


18 


93 





98 


78 


39 


1 


03 


-17 


89 





40 


(1 


1 


35) 


79 


TZ Boo 




74061 


15 


08 


09 


13 


+39 


58 


12 


86 


-89 


03 


1 


53 


58 
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1 


42 


6 


76 


1 


49 


-36 


70 


1 


51 


(1 


1 


28) 


80 
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74368 


15 


11 


5(1 


11 


+61 


51 


41 


26 
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96 


8 


04 


98 


81 


6 


44 


15 


32 


5 


18 


-60 


80 


1 


15 


(1 


1 


49) 


81 


BV Dra 


135421 


74370 


15 


11 


50 


36 


+61 


51 


25 


32 


-162 


57 


2 


98 


99 


13 


2 


47 


14 


86 


2 


56 


-61 


45 


1 


45 


(1 


1 


49) 


82 


FI Boo 


234224 


75203 


15 


22 


05 


97 


+51 


10 


55 


33 


-36 


87 


2 


40 


47 


43 


2 


29 


9 


52 


2 


10 


-30 


5 5 





72 


(1 


1 


35) 


83 


OU Ser 


136924 


75269 


15 


22 


43 


47 


+ 16 


15 


40 


74 


-387 


51 





91 


2 


80 





77 


17 


31 





95 


-04 


08 





41 


(1 


1 


18) 


84 


VZ Lib 




76050 


15 


31 


51 


76 


-15 


41 


00 


19 


-15 


61 


2 


55 


-1 


94 


1 


75 


4 


92 


1 


96 


-31 


11 


2 


30 


(1 


1 


35) 


85 
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76272 


15 


34 


45 


21 


+62 


16 


44 


28 


-255 


85 


1 


18 


15 


61 


1 


11 


6 


25 


1 


09 


-11 


38 


1 


14 


(1 


1 


18) 


86 
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141990 


77598 


15 


50 


32 


43 


+37 


50 


07 


56 


-74 


09 





92 


10 


93 





71 


11 


36 





85 


-4 


58 


1 


02 


(1 


1 


18) 


87 


AU Ser 






15 


56 


49 


47 


+22 


16 


01 


59 


-31 


40 


2 


50 


-48 


60 


2 


50 


5 


40 





68 


-62 


90 


1 


50 


(2 


5 


36) 


88 


V842 Her 






16 


06 


02 


34 


+50 


11 


12 


10 


-27 


90 


2 


50 


20 


70 


2 


50 


6 


08 





76 


-57 


98 


1 


55 


(2 


5 


26) 


89 


BX Dra 




78891 


16 


06 


17 


37 


+62 


45 


46 


10 


-7 


80 


1 


11 


7 


55 





92 


3 


38 


1 


15 


-26 


11 


3 


43 


(1 


1 


21) 


90 


V899 Her 


149684 


81191 


10 


35 


01 


96 


+33 


12 


47 


77 


39 


09 





72 


-39 


57 





55 


8 


06 





77 


-16 


84 


1 


30 


(1 


1 


35) 


91 
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150484 


81703 


16 


41 


20 


86 
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30 


27 


37 
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95 


1 


76 


20 


92 





86 


11 


84 


1 


17 


-42 


56 





85 


(1 


1 


21) 


92 


V918 Her 


151701 


82253 


16 


48 


24 


18 


+ 17 


07 


59 


63 


-15 


88 





88 


14 


27 





62 


8 


70 





93 


-25 


72 





74 


(1 


1 


21) 


93 


V921 Her 


152172 


82344 


16 


49 


31 


24 


+47 


06 


28 


81 


7 


40 





98 


5 
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91 


2 


39 





97 
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04 


1 


05 


(i;i;8) 


94 


V829 Her 






10 


55 


47 


88 
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10 


57 


60 


-11 


80 


3 


10 


-15 


60 


2 


60 


13 


53 





54 


-13 


40 


1 


00 


(4;f 


9 


11) 


95 


V2357 Oph 




82967 


16 


57 


16 


76 


+ 10 


59 


51 


38 





51 


1 


80 


-103 


60 


1 


11 


5 


37 


2 


01 


-19 


12 


1 


64 


(l;l;8) 


96 


AK Her 


155937 


84293 


17 


13 


57 


82 


+ 16 


21 


00 


61 


26 


24 


2 


23 


-47 


63 


1 


99 


10 


47 


2 


77 


-13 


00 


15.0 


(1 


1 


50) 


97 


V728 Her 






17 


18 


04 


46 


+41 


50 


40 


40 


-14 


40 


2 


50 


-10 


10 


2 


50 


2 


41 





30 


34 


40 





80 


(2 


5 


51) 


98 


GM Dra 


238677 


84837 


17 


20 


21 


88 


+57 


58 


27 


00 


77 


21 


1 


13 


-59 


05 


1 


02 


10 


16 





88 


9 


12 


1 


63 


(1 


1 


17) 


99 


V2377 Oph 


159356 


85944 


17 


33 


56 


05 


+08 


09 


57 


81 


-12 


87 


1 


00 


-59 


16 





72 


10 


09 


1 


22 


-25 


79 





38 


(1 


1 


35) 


100 


V535 Ara 


159441 


86306 


17 


38 


05 


55 


— 56 


49 


17 


29 


-76 


64 





90 


-49 


51 





50 


8 


87 





90 


-17 


60 





60 


(1 


1 


52) 


101 


V2388 Oph 


163151 


87655 


17 


54 


14 


17 


+ 11 


07 


50 


01 


-70 


81 





87 


-163 


11 





50 


14 


72 





81 


-25 


88 





52 


(1 


1 


17) 


102 


V566 Oph 


163611 
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17 


56 


52 


41 


+04 


59 


15 


32 
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04 


1 


17 


73 


92 





77 


13 
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1 


11 


-42 


95 


1 


13 
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66) 


103 
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164078 


87958 


17 


58 


05 


00 


+32 


38 


53 


02 


-34 


42 





51 


-28 


43 





46 


16 


25 





61 


3 


98 





46 


(1 


1 


17) 


104 


V508 Oph 




88028 


17 


58 


48 


62 


+ 13 


29 


46 


27 


-14 


62 


1 


88 
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07 


1 


38 


7 


68 


2 


14 


-36 


40 


2 


05 


(1 


1 


54) 


105 


EF Dra 






18 


05 


30 


48 


+69 


45 


15 


71 


-13 


60 


2 


70 


-24 


40 


2 


40 


5 


74 





72 


-42 


20 


3 


50 


(4 


5 


11) 


106 


V839 Oph 


166231 


88946 


18 


09 


21 


27 


+09 


09 


03 


63 


-32 


48 


1 


33 


3 


23 


1 


03 


8 


09 


1 


44 


-63 


76 





32 


(1 


1 


55) 


107 


cps Cra 


175813 


93174 


18 


58 


43 


38 


— 37 


06 


26 


48 


-132 


25 


1 


46 


-110 


45 





55 


33 


43 





92 


57 


90 


1 


20 


(1 


1 


56) 


108 


V401 Cyg 




95816 


19 


29 


20 


27 


+30 


24 


28 


48 


-4 


96 


8 


70 


-12 


59 


5 


99 


11 


55 


8 


04 


25 


53 


2 


14 


(1 


1 


17) 


109 


V2082 Cyg 


183752 


95833 


19 


29 


30 


07 


+36 


17 


14 


93 


20 


17 





53 


94 


83 





50 


11 


04 





56 


-34 


12 





58 


(1 


1 


21) 


110 


V417 Aql 




96349 


19 


35 


24 


12 


+05 


50 


17 


66 


-1 


90 


2 


86 


-24 


45 


1 


50 


7 


65 


2 


35 


-14 


20 


1 


40 


(1 


1 


11) 


111 


OO Aql 


187183 




19 


48 


12 


65 


+09 


18 


32 


38 


65 


70 


1 


00 


-7 


90 


1 


60 


6 


52 





82 


-45 


20 





90 


(4 


5 


57) 


112 


VW Cop 


197433 


101750 


20 


37 


21 


54 


+75 


36 


01 


46 


309 


23 


1 


12 


539 


47 





97 


36 


16 





97 


-7 


90 


1 


00 


(1 


1 


59) 


113 


LS Del 


199497 




20 


57 


10 


29 


+ 19 


38 


59 


23 


145 


40 


1 


90 


83 


40 


1 


80 


12 


71 


1 


59 


-25 


90 


1 


40 


(4 


5 


11) 


114 


V2150 Cyg 


202924 


105162 


21 


18 


10 


88 


+30 


35 


21 


58 


2 


58 


1 


27 


-6 


25 





91 


4 


73 


1 


59 


-12 


82 





45 


(1 


1 


35) 


115 


HV Aqr 






21 


21 


24 


81 


-03 


09 


36 


88 


-63 


50 


1 


90 


-92 


30 


1 


80 


6 


99 





87 


-3 


77 


1 


25 


(4 


5 


18) 


116 


V1073 Cyg 


204038 


105739 


21 


25 


00 


36 


+33 


41 


14 


94 


-4 


70 


2 


50 


-21 


00 


2 


50 


5 


44 





95 


-0 


80 


1 


10 


(2 


1 


60) 


117 


KP Peg 


204215 


105882 


21 


26 


41 


13 


+ 13 


41 


17 


75 


-2 


31 


2 


40 


-10 


65 


1 


20 


4 


37 


1 


67 


5 


52 


1 


03 


(1 


1 


21) 


118 


DK Cyg 




106574 


21 


35 


02 


67 


+34 


35 


45 


40 


-28 


37 


1 


24 


-19 


29 





93 


4 


42 


1 


78 


-5 


53 


1 


48 


(1 


1 


26) 


119 


BX Peg 






21 


38 


49 


39 


+26 


41 


34 


20 


7 


10 


2 


50 


-11 


90 


2 


50 


6 


13 





77 


25 


00 


3 


00 


(2 


5 


61) 


120 


V445 Cop 


210431 


109191 


22 


07 


10 


89 


+72 


22 


22 


26 


83 


91 





51 


43 


90 





42 


8 


95 





50 


40 


69 





95 


(1 


1 


21) 


121 


BB Peg 




110493 


22 


22 


56 


89 


+ 16 


19 


27 


84 


19 


28 


2 


66 


-24 


66 


1 


85 


3 


02 


2 


26 


-28 


90 


1 


40 


(1 


1 


11) 


122 


V335 Peg 


216417 


112960 


22 


52 


36 


18 


+10 


13 


13 


19 


166 


57 





86 


-89 


35 





65 


16 


26 





86 


-15 


41 





43 


(1 


1 


21) 


123 


SW Lac 


216598 


113052 


22 


53 


41 


66 


+37 


56 


18 


61 


88 


16 


1 


25 


10 


79 





82 


12 


30 


1 


26 


-14 


46 





90 


(1 


1 


62) 


124 


AB And 




114508 


23 


11 


32 


09 


+36 


53 


35 


11 


109 


02 


1 


55 


-53 


56 





94 


8 


34 


1 


48 


-27 


53 





67 


(1 


1 


21) 


125 
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220659 


115627 


23 


25 


25 


19 


+ 15 


41 


19 


14 


17 


43 





90 
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68 


7 
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92 


-8 


08 





89 


(1 


1 


14) 
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115819 


23 


27 


48 


39 


+04 


51 


23 


96 
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01 


2 


46 
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24 


1 


66 


16 


77 


2 


07 


-4 


30 


1 


80 


(1 


1 


63) 


127 


BD +14 5016 






23 


36 




37 


+ 15 


48 


06 


43 


19 


50 


1 


90 


2 


60 


1 


80 


3 


06 





38 


4 


20 





80 


(4 


5 


64) 
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EL Aqr 




117317 


23 


47 


18 


35 


-08 


05 


12 


09 


-7 


23 


2 


53 


1 


00 


1 


65 


4 


71 


1 


89 


12 


51 


2 


08 


(1 


1 


14) 
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U Peg 






23 


57 


58 


48 


+ 15 


57 


10 


08 


-46 


50 


1 


19 


-53 


88 





64 


7 


18 


1 


43 


-28 


50 


1 


60 


(1 


1 


65) 
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preferred to adopt an average uncertainty 2.5 mas yr -1 for 
them, which were determined by Hog et al. (2000) for the 
Hipparcos stars in general. On the other hand, Rucinski 
(2004) estimated that the absolute magnitudes computed 
by the period-brightness relation of W UMa systems are 
within the accuracy of ±0.25 mag. Therefore, a correspond- 
ing 12.5 per cent uncertainty for the photometric parallaxes 
were taken. 

Recent investigations, both visual and spectroscopic, 
indicated the presence of additional components for some 
W UMa binaries. Trigonometric parallaxes are frequently 
wrong for those systems as for the brighter systems with 
large uncertainties in the parallaxes. We have not pre- 
ferred to introduce additional uncertainties into the pub- 
lished data. However, for the reader is who want to estimate 
the bias introduced by such multiple system, the W UMa's 
with additional components are marked in the last column 
of Table 5 and are shown by empty symbols in the figures. 



2.2 Radial velocities 

All W UMa binaries have circular orbits. Thus the center of 
mass velocity can be obtained simply by fitting a sine curve 
to the measured radial velocities. However, except for the 
four systems (TW Cet, RZ Tau, TY Pup, AK Her), all W 
UMa stars in our list had well determined orbital parame- 
ters so that the center of mass velocities (7) were available 
with uncertainties. The center of mass velocities and the 
associated standard errors were listed in Table 1 with their 
sources. If there was more than one source available, the data 
in the latest source is taken since previous radial velocities 
were used in the later determined orbital parameters. Nev- 
ertheless, we have found the three systems (V523 Cas, AG 
Vir, V566 Oph) with multiple orbits which were determined 
independently from the independent data sets. For them, 
the weighted mean of the systemic velocities (7) and the 
weighted mean of the associated errors are adopted. Only 
the four systems (TW Cet, RZ Tau, TY Pup, AK Her) had 
measured center of mass velocities but uncertainties were 
not available. Their uncertainties were estimated from the 
scatter of component velocities on the radial velocity curve. 
All different types of errors have been transformed to stan- 
dard errors for the sake of consistency. 



3 GALACTIC SPACE AND VELOCITY 
DISTRIBUTIONS 

3.1 Sky and space distributions 

There exist no program similar to that of DDO, which con- 
tributed 54 per cent of radial velocities, in the southern 
hemisphere, where the number of observatories already rela- 
tively less so there is a strong bias in the sky distribution of 
the present W UMa systems. The celestial coordinates (a, 6) 
are plotted in Fig. la that the bias is obvious as there are 
103 systems (80 per cent) in the northern and 26 systems 
(20 per cent) in the southern hemisphere. Unfortunately, 
the bias is unavoidable because 1) The present work is es- 
tablished on the accumulated data in the literature and it 
is not vise to wait for another unknown period without a 
clue how southern systems will be increased. 2) Eliminating 
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Figure 1. Sky distributions of the present sample (a) at the celes- 
tial, (b) at the galactic coordinates. Circles: Hipparcos parallaxes; 
Triangles: the systems use Rucinski's calibration; Filled: normal 
W UMa; Empty: W UMa's with visible companions. 



some northern systems to establish isotropy is not feasible 
since it would greatly reduce the statistical significance of 
the sample. Due to the angle (63°) between the celestial 
equator and the galactic plane, distribution of the sample in 
the galactic coordinates (Fig. lb) appears to be less effected 
by the bias. Nevertheless, this is only an illusion because 
un-isotropy would not change by changing the coordinate 
system. 

In order to inspect the space distributions of the present 
W UMa sample, the Sun centered rectangular galactic coor- 
dinates (X towards Galactic center, Y Galactic rotation, Z 
north Galactic Pole) were calculated. The (X, Y, Z) space 
coordinates are given in Table 2. The projected positions on 
the Galactic plane (X, Y plane) and on the plane perpen- 
dicular to it (X, Z plane) are displayed in Fig. 2. The bias 
on the prime reference frames (X — Y , X — Z , and Y — Z) 
are not that obvious as Fig. la. Having 69 systems with 
X < and 60 system with X > 0, the Y — Z plane shows 
the least bias. Having 65 and 44 systems the northern and 
southern Galactic latitudes, Galactic plane (X — Y) shows 
the strongest bias (Fig. 2) but still it is not obvious as Fig. 
la. 

With a median distance of 137 pc, the sample contains 
relatively nearby systems and it can be considered within 
the Galactic thin disc. Occupying almost the same space in 
the Solar neighborhood, the CAB sample (Karatas, et al. 
2004), had a smaller median distance (98 pc) than W UMa 
sample. Thus, we were curious to compare space dispersions 
of the present W UMa sample and CAB sample. In order 
to avoid complications of far distant members, the space 
dispersions of both samples within 300 pc are calculated. 
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Table 2. Kinematic data of the W UMas. 



ID 


Name 


HD 


!(°) 


K°) 


X(pc) 


Y(pc) 


Z(pc) 


U (km 


a" 1 ) 


V (km 


s- 1 ) 


W (km 


s' 1 ) 


Code 


MG 


1 


AQ Tuc 


1372 


306 


.66 


-44 


97 


149 


-200 


-249 


-22 


76 


4.77 


-45.28 


4.28 


-9 


.93 


1.43 








2 


BH Cas 




119 


.08 


-3 


49 


-144 


259 


-18 


19. 


39 


3.76 


-15.02 


2.79 


-42 


.25 


6.37 








3 


DZ Psc 




118 


.30 


-41 


20 


-74 


137 


-136 


3 


.73 


2.42 


-3.05 


2.69 


25 


.60 


2.71 








4 


V523 Cas 




121 


.08 


-12 


.59 


-41 


67 


-18 


-43 


.78 


5.37 


-35.60 


4.21 


-28 


.31 


3.65 








5 


AQ Psc 


8152 


135 


.65 


-54 


.53 


-52 


51 


-101 


4 


.38 


0.70 


3.94 


1.57 


14 


.81 


0.81 








6 


AE Phe 


9528 


286 


.22 


-66 


.25 


5 


-19 


-45 


-20 


.12 


0.84 


-26.19 


1.53 


18 


91 


2.32 


1, 


1 


LA, IC 


7 


TW Cct 




190 


.59 


-75 


49 


-25 


-5 


-98 


-16 


.25 


7.52 


-20.14 


7.04 


-15 


.53 


24.24 


1> 1, 


2 


Cas, IC, LA 


8 


V776 Cas 




128 


.21 


7 


81 


-126 


160 


28 


-2 


16 


4.33 


-30.62 


3.90 


5 


.20 


2.95 








9 


QX And 




11 


.91 


11 


43 


365 


77 


75 


-24 


.12 


3.70 


-16.87 


4.41 


-17 


.36 


4.39 


1, 


1 


LA, IC 


10 


SS Ari 




143 


.57 


-35 


.94 


-138 


102 


-125 


-16 


.66 


2.87 


-24.10 


3.63 





.86 


2.18 


1, 


1 


IC, LA 


11 


GZ And 




137 


.72 


-15 


.87 


-116 


105 


-44 


5 


46 


1.90 


7.72 


1.89 


-4 


.07 


1.95 




1 


UMa 


12 


V376 And 


15922 


139 


.55 


-9 


.63 


-148 


126 


-33 


-51 


.49 


5.50 


-18.45 


5.92 


7 


.15 


2.03 








13 


EE Cot 


17613 


165 


.58 


-43 


.91 


-96 


25 


-95 


-19 


.66 


3.26 


-19.71 


4.03 


11 


.64 


2.87 








14 


UX Eri 




187 


.70 


-51 


.44 


-94 


-13 


-119 


-15 


.59 


3.79 


-14.45 


6.13 


-2 


.20 


3.67 


1, 


1 


IC, LA 


15 


EQ Tau 




168 


.13 


-24 


65 


-128 


27 


-60 


-81 


.50 


2.42 


-31.50 


5.71 


-14 


.52 


2.48 








16 


BL Eri 




205 


.10 


-40 


70 


-248 


-116 


-235 


-32 


.39 


3.65 


-23.90 


4.63 


-11 


.85 


4.04 








17 


YY Eri 


26609 


203 


.50 


-40 


04 


-39 


-17 


-36 


40. 


04 


2.22 


2.37 


0.57 


-20 


.96 


2.28 








18 


AO Cam 




152 


.16 


2 


.94 


-96 


51 


6 


3 


.61 


1.25 


-18.18 


1.92 


-3 


.46 


1.33 








19 


RZ Tau 


285892 


179 


21 


-18 


.72 


-165 


2 


-56 


1 


.85 


47.40 


9.78 


3.31 


-20 


.82 


17.25 








20 


DN Cam 


29213 


137 


.97 


17. 


30 


-158 


142 


66 


-12 


.77 


1.22 


9.48 


1.34 


-20 


.48 


4.48 




1 


UMa 


21 


V410 Aur 


280332 


170 


.06 


-4 


.65 


-206 


36 


-17 


-39 


.23 


3.88 


-2.64 


11.36 





.87 


6.13 








22 


V402 Aur 


282719 


172 


.79 


-6 


.45 


-141 


18 


-16 


-41 


.93 


0.96 


-1.58 


1.56 


-1 


.78 


1.09 








23 


AP Dor 


33474 


267 


.96 


-36 


.37 


-6 


-173 


-127 


-82 


.31 


17.16 


-30.34 


5.03 


33 


.94 


7.58 








24 


V1363 Ori 


289949 


201 


.09 


-23 


.44 


-90 


-35 


-42 


-22 


.53 


2.85 


-34.53 


5.65 


-15 


88 


1.27 








25 


ER Ori 




209 


.20 


-26 


.13 


-126 


-70 


-71 


-16 


64 


2.97 


-32.93 


2.67 


-20 


.00 


1.87 








26 


RW Dor 


269320 


278 


77 


-33 


.62 


14 


-92 


-62 


-23 


.25 


7.63 


-85.23 


6.89 





64 


8.40 








27 


V781 Tau 


248087 


182 


.21 


-0 


15 


-81 


-3 





-23 


.71 


1.30 


-12.75 


1.37 


-43 


.84 


4.83 








28 


AH Aur 


256902 


185 


.15 


7 


.31 


-160 


-14 


21 


-29 


16 


1.60 


-16.32 


4.75 


10 


.01 


2.54 








29 


QW Gem 


264672 


186 


19 


12. 


82 


-237 


-26 


54 


6 


.50 


7.06 


-30.72 


41.87 


6 


.47 


12.26 








30 


V753 Mon 


54975 


218 


.61 


2 


.56 


-149 


-119 


9 


-26 


.33 


0.78 


-25.43 


0.81 


-12 


.98 


3.04 




2 


IC 


31 


TY Pup 


60265 


236 


.04 


-0 


.72 


-120 


-178 


-3 


-51 


62 


17.50 


4.56 


25.11 


-1 


.09 


1.11 








32 


FG Hya 




221 


.03 


22 


83 


-238 


-207 


133 


61 


63 


33.93 


-78.64 


53.82 


-37 


89 


23.38 








33 


TX Cnc 




206 


.62 


32. 


17 


-145 


-73 


102 


-42 


.09 


3.67 


-22.51 


3.60 


-24 


.46 


5.27 








34 


AH Cnc 




215 


.68 


31. 


99 


-574 


-412 


441 


-6 


.11 


7.33 


-35.39 


9.38 


-28 


54 


9.59 








35 


UV Lyn 




184 


.74 


41. 


51 


-91 


-8 


81 


-27 


.58 


5.69 


12.21 


2.45 


-30 


.62 


6.15 








36 


FN Cam 


79886 


135 


15 


34 


44 


-162 


161 


156 


-30 


.57 


5.05 


-22.65 


8.09 


19 


53 


3.45 




1 


LA 


37 


EZ Hya 




245 


.81 


25 


69 


-68 


-151 


80 


-20 


97 


3.38 


-14.28 


1.31 


-11 


.86 


3.29 


1, 


2 


LA, IC 


38 


S Ant 


82610 


258 


.55 


Hi. 


61 


-14 


-71 


21 


-30 


51 


1.82 


8.15 


1.70 


-12 


.40 


0.80 








39 


W UMa 


83950 


158 


.92 


45. 


90 


-32 


12 


36 


21. 


01 


0.82 


-13.19 


0.45 


-17 


.29 


0.91 








40 


AA UMa 




173 


.25 


49. 


21 


-203 


24 


237 


20. 


99 


3.09 


-9.41 


3.81 


-26 


.42 


2.82 








41 


R.T Lmi 




190 


.74 


50. 


04 


-203 


-38 


252 


13. 


20 


3.29 


-6.83 


3.98 


-4 


.65 


2.82 








42 


XY Leo 




217 


80 


49. 


74 


-32 


-25 


48 


45. 


60 


2.46 


6.51 


1.72 


-32 


.13 


1.70 








43 


XZ Leo 




218 


.44 


49. 


80 


-195 


-155 


295 


-19 


18 


3.99 


0.06 


1.54 


-18 


.35 


3.03 


1, 


1 


LA, Cas 


44 


Y Sex 


87079 


238 


.49 


41. 


94 


-43 


-70 


74 


-0 


.95 


1.88 


-7.13 


1.58 


6 


.14 


1.87 








45 


ET Leo 


91386 


222 


.34 


56 


81 


-29 


-27 


60 


-17 


.73 


1.15 


-13.03 


0.67 


10 


.95 


1.01 








46 


UZ Leo 




230 


.28 


56. 


60 


-56 


-68 


133 


-8 


.66 


3.55 


-1.22 


1.46 


-14 


.91 


2.32 




1 


UMa 


47 


EX Leo 


93077 


226 


.62 


58. 


85 


-36 


-38 


87 


1 


.33 


0.77 


-15.00 


2.24 


-19 


.42 


1.49 








48 


VY Sex 


93917 


253 


.81 


48. 


30 


-25 


-87 


102 


15. 


60 


2.67 


-56.00 


5.65 


-21 


.28 


4.39 








49 


AM Leo 




241 


50 


59. 


02 


-19 


-35 


66 


1 


26 


1.62 


-12.59 


3.08 


-0 


05 


1.87 




1 


Cas 


50 


VW Lmi 


95660 


198 


.56 


66. 


04 


-48 


-16 


114 


5 


48 


0.98 


-0.49 


0.42 


7 


.51 


0.80 


1, 


2 


Cas, UMa 


51 


AP Leo 




249 


.13 


56. 


47 


-24 


-62 


101 


67. 


11 


12.50 


2.03 


2.59 


-13 


.45 


2.39 








52 


HN UMa 




179 


28 


67. 


35 


-66 


1 


159 


63. 


15 


11.73 


-11.70 


2.90 


-13 


80 


4.94 








53 


AW UMa 


99946 


199 


.11 


71. 


93 


-19 


-7 


63 


7 


01 


0.54 


-65.57 


4.02 


-22 


.78 


1.55 








54 


TV Mus 


310730 


295 


.42 


-2 


.99 


99 


-209 


-12 


34. 


61 


6.06 


10.03 


2.97 





.76 


3.72 








55 


V752 Ccn 


101799 


287 


.63 


25. 


01 


29 


-91 


44 


-8 


12 


2.97 


-39.68 


2.37 


-5 


.22 


2.77 








56 


AG Vir 


104350 


260 


.60 





72 


-33 


-197 


3 


11. 


70 


1.43 


-13.09 


1.99 


-10 


.20 


0.94 








57 


HX UMa 


104425 


152 


.16 


71. 


14 


-43 


23 


142 


51. 


86 


21.15 


4.70 


3.77 


-5 


.90 


5.94 








58 


CC Com 




238 


.04 


79. 


97 


-7 


-12 


79 


-41 


.89 


5.34 


-13.52 


2.07 


-16 


.34 


1.51 








59 


AH Vir 


106400 


271 


.52 


72. 


39 


1 


-28 


88 


42 


13 


11.89 


-30.79 


8.29 


-3 


.19 


2.86 








60 


II UMa 


109247 


128 


.65 


62. 


13 


-58 


72 


175 


-29 


.02 


10.65 


-32.09 


10.61 


-4 


.69 


1.45 


1, 


1 


Hya, IC 


61 


RW Com 




217 


.61 


85. 


87 


-5 


-4 


87 


-33 


.71 


7.99 


-37.02 


8.48 


-56 


.70 


1.41 








62 


RZ Com 




257 


.75 


84. 


70 


-4 


-19 


212 


21. 


35 


3.06 


-0.24 


1.39 


-1 


.42 


1.67 








63 


SX Crv 


110139 


299 


.25 


43. 


99 


32 


-58 


64 


19. 


62 


1.76 


1.69 


1.06 


5 


.01 


0.72 








64 


BI CVn 




108 


.86 


80. 


19 


-10 


30 


186 


3 


.13 


1.55 


-16.11 


4.86 


-2 


.52 


1.97 




1 


Cas 


65 


KZ Vir 


114726 


315 


.40 


65 


02 


46 


-45 


139 


-5 


49 


0.99 


-1.82 


0.67 


-3 


.29 


0.48 


1, 


1 


UMa, Cas 



M 44 
M 67 



Wolf 630 



Indeed, with the space dispersions of 5X = ±78, 5Y = ±90, 
and 5Z = ±87 pc, the W UMa sample seems a little more 
dispersed than the CAB sample with the space dispersions 
of SX = ±77, SY = ±75, and SZ = ±73 pc. 

Considering the fact that CAB's being on average 
brighter than the less common W UMa's, it has been ex- 
pected the opposite, since CAB's would be visible at longer 
distances. Obviously, the observational selection is more af- 
fected by the easy recognition of W UMa's rather than the 



higher absolute brightness of CABs. Apparently, the less 
bright contact binaries were identified and studied more 
than the CAB's in similar distances. 

3.2 Galactic space velocities 

Galactic space velocity components (U, V , W) were com- 
puted by applying the algorithm and the transformation 
matrices of Johnson & Soderblom (1987) to the basic data: 
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Table 2 - continued 



ID 


Name 


HD 




b(°) 


X(pc) 


Y(pc) 


Z(pc) 


U (km 


s- 1 ) 


V (km 


s- 1 


) 


W (km 


s- 1 


) 


Code 


MG 


OC 


66 


KR Com 


115955 


339 


81 


78.52 


14 


-5 


75 


13 


.68 


1.03 


-24.99 


1 


.71 


-12 


.35 





49 








67 


HT Vir 


119931 


335 


.97 


64.42 


26 


-11 


59 


-32 


.25 


4.21 


-19.12 


4 


.16 


-15 


.40 


1 


.21 








68 


XY Boo 




8 


.30 


75.08 


87 


13 


329 


-54 


.47 


31.71 


10.42 


6 


.48 


16 


05 


8 


36 








69 


V757 Con 


120734 


316 


25 


24.68 


46 


-44 


29 


6 


.12 


2.73 


-66.49 


3 


15 


12 


.13 


1 


28 








70 


RR Cen 


124689 


314 


13 


3.15 


71 


-73 


6 


-28 


.95 


1.77 


-8.09 


1 


.73 


-2 


.83 





34 








71 


VW Boo 




1 


29 


65.15 


122 


3 


264 


54 


78 


24.84 


-46.53 


25 


.75 


1 


.43 


11 


26 








72 


NN Vir 


125488 


351 


.55 


60.16 


52 


-8 


91 


-16 


.14 


1.68 


-24.65 


3 


.03 





.02 





89 


1 


Cas 




73 


EF Boo 


234150 


90 


.81 


59.53 


-1 


85 


144 


2 


.98 


1.23 


-40.40 


6 


.03 


8 


.54 


3 


.66 








74 


CK Boo 


128141 


1 


.46 


59.62 


79 


2 


135 


89 


.49 


14.89 


-13.67 


3 


02 


-9 


.34 


8 


.71 








75 


GR Vir 


129903 


346 


.23 


46.34 


36 


-9 


38 


-65 


.96 


1.31 


6.30 





.48 


-36 


.32 


1. 


.19 








76 


AC Boo 




79 


26 


58.76 


13 


67 


113 


-10 


53 


1.55 


-5.11 


1 


.32 


-6 


.05 


1 
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celestial coordinates (a, 5), proper motion components (/j, a , 
Us), radial velocity (7) and the parallax (tt) of each star in 
Table 1, where the epoch of J2000 coordinates were adopted 
as described in the International Celestial Reference System 
(ICRS) of the Hipparcos and the Tycho Catalogues. The 
transformation matrices use the notation of the right handed 
system. Therefore, the U, V, W are the components of a ve- 
locity vector of a star with respect to the Sun, where the U 
is directed toward the Galactic center (I = 0°, b = 0°); the V 



is in the direction of the galactic rotation (I = 90°, b = 0°); 
and the W is towards the north Galactic pole (b = 90°). 

Uncertainties of (U, V, W) space velocity components 
have been computed by propagating the uncertainties of 
the input data (proper motions, parallax and radial veloc- 
ity) by the algorithm also by Johnson & Soderblom (1987). 
The arithmetic mean of the uncertainties for (U , V, W) 
are: SU = ±5.52, SV = ±4.67 and SW = ±3.75 km s _1 . 
By inspecting the individual errors of space velocity vec- 
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Figure 2. The space distributions with respect to the Sun. X, Y, 
and Z are heliocentric galactic coordinates directed towards the 
Galactic Center, Galactic rotation and the North Galactic Pole. 
Symbols mean the same as Fig. 1. 



Figure 3. (a) Velocity dispersions on the U — V plane, (b) propa- 
gated errors of some selected systems. Symbols in (a) are the same 
as Fig. 1 and symbols in (b) are self indicated. The velocities are 
heliocentric. The position of LSR is marked by +. 



tors (S = VU 2 + V 2 + W 2 ), we have found only 11 (8.5 per 
cent) systems with the errors bigger than ±20 km s _1 . If 
those systems were removed from the sample, the averages 
would reduce to SU = ±4.04, SV = ±3.12, and SW = ±2.68 
km s~ x . Thus, it can be claimed confidently that most of the 
sample stars have the velocity uncertainties smaller than the 
velocity dispersions calculated. However, in comparison with 
the CAB sample (Karatas, et al. 2004), the average uncer- 
tainties of the present W UMa sample appear nearly two 
times bigger. Obviously, this is the result of the lower accu- 
racy of W UMa radial velocities because of shorter orbital 
periods and higher rotation rates. 

Having relatively less accurate space velocities, we were 
curious to check the effect of the differential Galactic rota- 
tion. Thus, the first-order Galactic differential rotation con- 
tributions to the U and V components were computed as 
described in Mihalas & Binney (1981). The W velocities are 
not affected in this first order approximation. We have found 
27 stars (21 per cent) with the contribution from the Galac- 
tic differential rotation being greater than the uncertainty 



of the U component of the space velocity. The effect on the 
V component is rather small that the effect of Galactic dif- 
ferential rotation was found smaller than the uncertainties 
in V . Nevertheless, the Galactic differential rotation correc- 
tion was applied to all of stars in the present sample. The 
corrected U, V, W are given in Table 2, together with the 
propagated standard errors. 

There is no direct way to estimate how the un-isotropic 
sky distribution (bias shown in Fig. la) effects the W UMa 
space velocities. However, since present sample contains only 
nearby systems within the Galactic thin disc and since first 
order Galactic differential rotation correction were applied, 
the bias seems to be minimized perhaps disappeared in the 
space velocity distribution. This is because the velocity dis- 
tribution present W UMa sample has a great resemblance 
to the velocity distribution of the CAB sample (Karatas, et 
al. 2004) which were considered homogeneously and isotrop- 
ically distributed in the solar neighbourhood. Distribution 
of the space velocity components on the U — V plane is dis- 
played in Fig. 3a. 
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In order to display reliability of computed U, V , W 
velocities, Fig. 3b shows error bars some selected systems as 
1) five system which are the most distant, 2) five systems 
with the biggest space velocity vector (S), 3) five systems 
which are the most dispersed on the U — V diagram. It is 
also clear on Fig. 3b that propagated uncertainties, even the 
extreme cases as above, appear smaller than the dispersion 
thus the space velocities (U, V , W) and their dispersion can 
be considered realible. Nevertheless, in order to minimize the 
contribution of the space velocity errors, the dispersion were 
calculated with respect to the LSR rather than arithmetic 
mean of the space velocities. The position of the LSR is 
obtained by subtracting the Sun's velocity (U,V,W)q = 
(9, 12, 7) km s _1 by Mihalas & Binney (1981) from the frame 
of (U, V, W) velocities. 



3.3 Population analysis 

As in the CAB sample (Karata§ et al. 2004), we have used 
the parameter f K = (l/300)(1.0u 2 + 2.5i> 2 + 3.5to 2 ) 1/2 to 
determine the possible metal poor binaries kinematically as 
suggested by Grenon (1987) and Bartkevicius, Lazauskaite, 
& Bartasiute (1999). Here, the u, v, w velocities represent 
space velocity components with respect to the LSR. Statis- 
tically, the stars fx < 0.35 belong to the thin disc, the stars 
with 0.35 < < 1-00 belong to the thick disc. The stars 
with fa > 1 belong to the halo population. 

According to f K analysis, the vast majority (93 per 
cent) of our sample are thin disc stars. The possible thick 
disc members are found to be 7 per cent in the current sam- 
ple. The thick disc candidates are FU Dra, AP Dor, RW 
Dor, FG Hya, OU Ser, AU Ser, V921 Her, V2357 Oph, 
and VZ Psc. Spectroscopic metal abundances of W UMa 
stars appear unachievable due to fast rotation and Mend- 
ings in their spectra. Photometric indices are also entirely 
unreliable due to the elevated magnetic activity of the W 
UMa systems. Nevertheless, photometric metal abundance 
of BW Dra has been determined by Marsakov & Shevelev 
(1995) ([m/H] — —0.56 dex) and recently by Nordstrom 
ct al. (2004) ([m/H] — —0.69 dex) from the observations 
of Stromgren photometry. A 10.33 Gyr age is estimated for 
it by Marsakov & Shevelev (1995) according to the revised 
Yale isochrones ages (Green, Demarque, & King 1987). If 
those metal abundances and isochrone ages are correct, with 
the /k kinematic value of 0.33, BW Dra appears as a star 
which belongs to the thick disc population although it's Jk 
value is little less than the lower limit. 

Karata§ et al. (2004) had found 7 per cent of CAB bina- 
ries are possible members of the thick disc population. It is 
indeed interesting that the same percentage of thick disc W 
UMa binaries confirmed in this study coincides with nearly 
the same ratio, 6 per cent of solar neighborhood stars in gen- 
eral belongs to the thick disc population which were found 
by Buser, Rong, & Karaali (1999), Siegel et al. (2002), and 
Karaali, Bilir, & Hamzaoglu (2004). 



4 DISCUSSION 

4.1 Kinematics of contact binaries 

Having only the main sequence components from the spec- 
tral types A to K, the present sample of W UMa binaries 
appear less heterogeneous than the CAB binaries with the 
components F to M spectral types of all combinations of 
the luminosity classes from V to II. Having orbital periods 
from 1 days to several hundred days, the CAB sample is 
also considered heterogeneous (Karatas, et al. 2004) regard- 
ing to the period ranges since different orbital periods may 
represent different evolutionary paths (Plavec 1968, Thomas 
1977). But, with a very limited range of short orbital periods 
(0.22 < P < 1.13 days), the W UMa sample can be claimed 
homogeneous in a sense that all binaries with main sequence- 
components are in the contact stage. However, the velocity 
distribution of W UMa sample on the U — V plane (Fig. 
2) appears heterogeneous as the CAB sample, that kine- 
matically younger and older sub systems occupy the same 
velocity space. With the dispersions of 36.5, 26.2, 19.5 km 
s _1 at U, V, and W velocity components, it is possible to 
assign a 4.43 Gyr kinematical age to the whole sample ac- 
cording to the kinematical tables of Wielen (1982). Karata§ 
et al. (2004) assigned 3.86 Gyr age to the field CAB binaries. 
Thus, even before selecting out possible moving group (MG) 
members, which indeed seems to exist from the appearance 
of Fig. 3, the W UMa sample is little older than the CAB 
sample. 

Guinan & Bradstreet (1988) investigated whether cor- 
relations exist between the space velocities and the physical 
properties of the W UMa systems, such as orbital period 
(P), mass ratio (q), system type (A or W), filling factor (/), 
ultraviolet excess, and metallicity index. A weak correla- 
tion appeared possible in the sense that the space velocities 
tend to increase with decreasing orbital period. Also, a weak 
correlation appeared as binaries with large / (over-contact) 
tend to have higher space velocities. But, no correlation had 
been found between the space velocities and the mass ratio, 
the space velocities and the binary type (A or W), the space 
velocities and the metallicity index. 

In order to search kinematically different sub groups, 
here we have preferred to divide the whole sample into two 
sub samples according to a criterion, then to compare their 
U — V diagrams. The preliminary criteria and the U — V 
diagrams are summarized and displayed in Fig. 4. 

W UMa binaries with spectral types later than F7 tend 
to have larger dispersions (Fig. 4a and b) than W UMa's of 
earlier spectral types. This could be explained by the fact 
that the later spectral types on the H-R diagram contain 
older systems since evolution into the main-sequence and 
the duration on the main sequence is longer. However, it 
is indeed extraordinary that kinematically young and old 
systems (small and large dispersions) may well exist in both 
groups. Thus, this criterion fails to provide kinematically 
homogeneous sub samples. 

W UMa stars are classified into A- or W-type systems 
from their light curves and radial velocity curve (Binnendijk 
1970). The W-type systems are those like W UMa itself in 
which the hotter component (the star eclipsed at the pri- 
mary minimum) is the smaller and the less massive. Com- 
parison of the U — V diagram between the A- and W-type in 
Fig. 4 confirms the Guinan & Bradstreet (1988) that both 
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Figure 4. Comparison of the space velocity dispersions of various sub groups on the U — V diagram, (a) earlier spectral types A2 to F7; 
(b) later spectral types F8 to K5; (c) A type; (d) W type; (e) mass ratio q < 0.32; (f) mass ratio q > 0.32; (g) over-contact parameter 
/ < 0-2; (h) over-contact parameter / > 0.2. 
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Figure 5. The space velocities is plotted against the fillout factor 
(/) for W UMa binaries. 



groups are indistinguishable kinematically, that younger and 
older systems are equally likely to exist in both sub groups 
of A or W. 

If a few high dispersion binaries are removed from the 
small mass ratio systems, then both sub groups (Fig. 4e and 
f ) of mass ratio criterion will be indistinguishable kinemat- 
ically. Like the two criteria before, the mass ratio criterion 
also fails to distinguish younger and older systems homoge- 
neously. 



The degree of over-filling (/) of the inner critical equipo- 
tential surface is an important parameter which may indi- 
cate a stage of an evolution for a contact binary which would 
coalesce into a single star. FK Com is believed by some to 
be the product of binary star coalescence (Bopp & Rucinski 
1981, Bopp & Stencel 1981). The value of over-filling factor 
/, which depends on how it is defined, changes from the min- 
imum, when the system is just in contact, to the maximum 
when components fillout their outer critical envelopes. Dif- 
ferent authors define f differently, e.g. Guinan & Bradstreet 
(1988) define f as 1 < / < 2, but Maceroni & van't Veer 
(1996) define f as < / < 1. Using the same definition as 
Maceroni & van't Veer (1996), we have found 109 systems 
with (/). The systems with smaller / are compared to the 
systems with large / on the U — V diagram (Fig. 4g and h). 
The systems with larger / values seem to have larger disper- 
sions, so they appear kinematically older than the systems 
with smaller /. Although this is a natural consequence of the 
theory of coalescence into a single star, and it may appear 
to confirm Guinan & Bradstreet (1988), with a careful look 
at the U — V diagrams (Fig. 4g and h) one may notice the 
two sub groups are not really homogeneous. Kinematically 
older and younger systems appear to exist in both groups. 

In order to make a direct comparison to the Fig. 9 
of Guinan & Bradstreet (1988), the space velocities (s = 
\/u 2 + v 2 + w 2 ) in our sample is plotted against the / values 
in Fig. 5. According to this presentation, the high velocity 
stars may well exist at low or high values of /. Thus it can 
be concluded that / criterion, like the criteria before (spec- 
tral type, binary type (A or W), mass ratio) is not a good 
criterion to differentiate kinematically older and younger W 
UMa binaries. 

This result, however, is consistent with the very brief 
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lifetime of a contact stage, which was predicted by Guinan 
& Bradstreet (1988) as 0.1 < t cont act < 1 Gyr, and the 
kinematical ages of the sample as a whole which is 4.43 
Gyr. The dispersions on the U — V diagram are produced 
in a time scale not only covering the contact stage but also 
covering the pre-contact stages of all. Obviously the contact 
stage for many W UMa stars is too short to be effective on 
the observed dispersions. Therefore, it is possible that a W 
UMa system with a longer pre-contact evolution may reach 
to a contact configuration at kinematically larger ages that 
the / value can be very small. But, on the other hand it is 
also possible for a W UMa system, with a much shorter pre- 
contact stage (small kinematical ages), to have larger value 
for the /. Therefore / value criterion fails to discriminate 
between kinematically young and old systems. 

The present W UMa sample were said above to be ho- 
mogeneous in a narrow range (0.22 < P < 1.13 days) of the 
short orbital periods. However, the correlations of orbital 
periods with the dispersions on the U — V diagram is well 
presented in Fig. 6 in a sense that the short period W UMa 
systems have larger dispersions. The sub groups according 
to the total masses (M = Mi + M2) display a similar corre- 
lation that less massive systems have larger dispersions than 
more massive systems. 

All sub groups discussed above are summarized numer- 
ically in Table 3. The numerical values (velocity averages 
and dispersions) of the sub groups confirm eye inspection of 
the dispersions on U — V diagrams on the Fig. 4 and Fig. 
6. Indeed the grouping according to spectral type, systems 
type (A or W), mass ratio, and / value in Fig. 4 do not 
show very much differences in the dispersions. But group- 
ing according to orbital period and total mass (Fig. 6) is 
confirmed strongly in Table 3. Using the tables of Wielen 
(1982), a kinematical age was assigned to each hypothetical 
sub-group. The ages represent a mean kinematical age of 
all the stars in a group. The age to be meaningful depends 
on the meaningfulness of the group. The sub groups using 
orbital period and total mass criteria are indeed meaningful 
because the decrease of orbital period and total mass, which 
is established for CAB (Karata§ et al. 2004), appears to be 
a process continuing also in the contact stage. Because of 
a relatively very short contact stage, there is no guarantee 
about the sub groups of total mass and orbital period being 
homogeneous. One way of selecting out younger population 
W UMa binaries among the sub samples according to or- 
bital period, is to pick out possible moving group members 
as done in CAB sample by Karata§ et al. (2004). 

4.2 Possible MG members among W UMa 
systems 

Moving groups (MGs) are kinematically coherent groups of 
stars that share a common origin. Eggen (1994) defined a 
supercluster of stars gravitationally unbound in the solar 
neighborhood, but sharing the same kinematics while occu- 
pying the extended regions in the Galaxy. Therefore, a MG, 
unlike the well known open clusters covering a limited re- 
gion, can be observed all over the sky. The kinematical crite- 
ria, which were originally defined by Eggen (1958a, b, 1989, 
1995), for determining the possible members of the best- 
documented MGs are summarized by Montes et al. (2001a, 
b). Evidence has been found that many young and active 



late-type binaries can be the members of some young mov- 
ing groups (Jeffries 1995, Montes et al. 2001b, King et al. 
2003). Indeed, possible moving group members determined 
by Eggen's kinematical criteria have been proved to be very 
useful for separating kinematically heterogeneous CAB sam- 
ple into two sub samples representing younger and older 
populations (Karatas, et al. 2004) better than the classical 
approach of pre determined sub groups with a dividing line 
of a chosen parameter. 

The difficulty of separating kinematically young and old 
populations in the velocity space is obvious. The dispersions 
increase with age but there are always some stars naturally 
occupying the regions near the LSR. It is, therefore, not 
safe to pick stars randomly near the LSR and then to form 
a kinematically young group with them. Eggen's kinemat- 
ical criteria, however, objectively select some stars which 
have space velocity vector parallel to the converging point 
of each pre-determined and well known MGs within an ac- 
ceptable uncertainty. There are two criteria that one sets a 
range of the direction of the test star's proper motion and 
other puts limits on the test stars radial velocity vector. 
Fulfilling one of the criteria makes the test star a possi- 
ble member. But, fulfilling both criteria does not guarantee 
membership since there is always a possibility that the same 
velocity space could be occupied by the MG members and 
the non-members. Further independent criteria implying a 
common origin and same age as the member stars are needed 
to confirm the true membership. Therefore, we should al- 
ways remember that Eggen's criteria determines only the 
possible members of a MG. The parameters of the five best- 
documented MGs and the possible membership criteria to 
them have been summarized by Karatas, et al. (2004) for 
the CAB binaries. Here we apply the same criteria to the 
stars in the present W UMa sample and 28 systems were 
found to satisfy at least one of the criteria for one of the 
five MGs which were listed in Table 4. The possible moving 
group members among the W UMa binaries in the present 
sample is marked in the latest columns of Table 2 with the 
names of MG associated. 

After determining the possible moving groups, the 
whole sample has been divided into three sub groups. The 
first one contains the possible moving group members and 
named MG. The second group contains only five stars which 
are previously known to be the members of well known open 
clusters. The third one is called 'Field Contact Binaries' 
(FCB) which contains field W UMa systems free from the 
possible MG and the known open cluster members in the 
solar neighborhood. Fig. 7 compares the U — V diagram of 
the MG and the FCB. The detailed kinematics and implied 
kinematical ages of MG and FCB are included in Table 3. 

The velocity dispersions of field contact binaries (FCB) 
imply a mean kinematical age of 5.47 Gyr. The mean kine- 
matical age is found to be 500 Myr for the possible MG 
members. Although, this age appears to be consistent with 
the ages of MG in Table 4, it may not represent the true 
mean kinematical ages of all possible MG members. This 
is because, the ages in Table 3 are estimated from the dis- 
persions with respect to the LSR. In reality, each MG has 
its own center of dispersion which is slightly different than 
the LSR. Therefore, we should keep in mind that members 
should be at the age of each MG listed in Table 4. 
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Figure 6. Comparison of the space velocity dispersions of period and total mass sub groups on the U — V diagram, (a) 0.2 < P < 0.4 
days; (b) 0.4 < P < 0.5 days; (c) P > 0.5 days; (d) M to t < 1.6OM ; (e) 1.60 < M to t < 1.95M Q ; and (f) M tot > 1.95M Q . 



Table 3. Kincmatical data and the predicted kincmatical ages of the sub groups of W UMa. 



Data 
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< U 
(km s 


> 
-1 




< V 
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> 
-1 




< W 
(km s 


> 
-1 






"U 
km s" 


-1 


) 
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-1\ 


(km s 


-1\ 


at 
(km 


ot 

— 1\ 


Age 

(Gyr) 


All 
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-12.11 





62 


-21.19 





62 


-8.52 





41 


36 


50 


7 


06 


26 


23 


7.09 


19 


46 


4.65 


49.0 


11.0 


4.43 


A2-F7 Spectral typo 


64 


-12.28 





96 
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51 


-6.16 





55 


28 


94 
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68 


21 


34 


4.11 


17 


22 


4.43 


39.9 
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2.83 


F8-K5 Spectral type 
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12 
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60 


42 
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6 


42 


30 
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21 


46 
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12.1 


5.88 


A Type 
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88 
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39 


86 


8 
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24 
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15 
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71 
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21 


4 


80 


28 


41 
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22 


85 
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10.0 
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96 


-3.74 





49 


40 


97 


7 


29 


29 
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22 


08 
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80 
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67 


31 


15 


6 


80 


22 


52 
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16 


26 
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31 
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41 


5 


30 


24 


36 
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18 
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Table 4. The parameters of the best documented moving groups. 



Name 


Age 


(U, V, W) 


V T * 


C. P.** 




(Myr) 


(km s _1 ) 


(km s~ 1 ) 


(a h , 5°) 


Local Association 


20 - 150 


(-11.6,-21.0,-11.4) 


26.5 


(5.98,-35.15) 


(Pleiades, a Per, M34, 










<5 Lyr, NGC 2516, IC 2602) 










IC 2391 Supercluster 


35 - 55 


(-20.6,-15.7,-9.1) 


27.4 


(5.82,-12.44) 


(IC 2391) 










Castor MG 


200 


(-10.7,-8.0,-9.7) 


16.5 


(4.75,-18.44) 


Ursa Major Group 


300 


(14.9,1.0,-10.7) 


18.4 


(20.55,-38.10) 


(Sirius Supercluster) 










Hyades Supercluster 


600 


(-39.7,-17.7,-2.4) 


43.5 


(6.40,6.50) 


(Hyades, Praesepc) 











* Space velocity vector of MG as a whole (V T = VU 2 + V' 2 + W 2 ). 

** Direction of space velocity vector (C. P.: Converging Point in celestial coordinates). 
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Figure 7. Distribution of (a) possible MGs and (b) FCBs on the 
(U, V) diagram. The position of LSR is marked by +. Symbols 
are the same as Fig. 1. 



4.3 Comparing physical parameters between FCB 
and MG 

The kinematical criteria select out kinematically young W 
UMa systems without a dividing line according to any physi- 
cal parameter of W UMa characteristics. However, this does 
not mean all kinematically young W UMa's are removed 
from the main sample when forming the group MG. There 
could be systems as young as the MG group in the FCB. 
Those are the systems which fail to satisfy kinematical cri- 
teria although they are young. Therefore young and old all 
sorts of W UMa system may exits in FCB. Thus, the 5.47 
Gyr kinematical age represents an average age of all stars in 
FCB. On the other hand, One must remember that the MG 
group stars are not young because of their small dispersions 
on the velocity space. On the contrary, their dispersions are 
small because they are young. They are young primarily be- 
cause they are possible members of young moving groups 
in Table 4. Nevertheless, the word "possible" implies a lim- 
ited number of non-members may well exist in MG which 
we believe not changes the statistic of young W UMa phe- 
nomenon. 

Comparing the MG group and field CAB's, Karatas, et 
al. (2004) found some important clues indicating mass loss 
and orbital period decrease in the binary evolution. A simi- 
lar study is intended here to study the contact stage of the 
binary evolution. The physical parameters of W UMa bina- 
ries are listed in Table 5. The columns are self explanatory 
to indicate order number, name, apparent brightness and 
color (B — V), spectral type, binary type (A or W), orbital 
period and inclination, masses of components, mass ratio, 
and over-contact parameter /. The data was primarily col- 
lected from the same literature where the radial velocities 
were taken. Moreover, the orbital inclinations and / values 
are taken from Pribulla, Kreiner, & Tremko (2003) and Se- 
lam (2004). 

The distribution according to spectral types are com- 
pared between the MG (younger) group and FCB (older) 
group in Fig. 8. It is interesting that both groups have a 
similar shape distribution. A slightly increased percentage 
of A type stars in MG group appears to be the only dif- 
ference. Fig. 9 compares the orbital period distributions be- 
tween MG and FCB. If one of the stars with an orbital period 
greater than 1 day in MG group is not counted, both groups 
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cover nearly the same orbital period range (0.2 < P < 0.9 
days). Both groups can be said to have a peak at P = 0.4 
days. With 27 members the MG group has an average pe- 
riod < P >= 0.5258 days. With 97 members FCB group 
has an average period < P >= 0.4275 days. Comparing the 
shape of both histograms we can conclude that the relative 
number of longer period (P > 0.4 days) systems decreased 
and the relative number of shorter period systems (P < 0.4 
days) increased in the older population sample (FCB). 

These changes may appear as evidences of the period 
decrease in the contact binary evolution. Karata§ et al. 
(2004) assumed the period histogram of possible MG mem- 
bers as the initial period distribution of field CAB. But, 
here, the situation is totally different. If the contact stage 
time scale is indeed short as estimated to be within the range 
0.1 < t cont act < 1 Gyr by Guinan & Bradstreet (1988), MG 
contact binaries cannot be taken as an initial distribution of 
the FCB, because MG contact binaries will not be able to 
survive up to the mean kinematical age of FCB. 

The mean kinematical age of FCB is 5.47 Gyr (Table 
3). Karata§ et al. (2004) assigned 3.86 Gyr age for the field 
CAB which are potential progenitors of field contact binaries 
(Demircan 1999). Therefore, if there is a continuous forma- 
tion of CAB and evolution out of the CAB stage, some (or 
maybe most) of them become W UMa systems, and simi- 
larly, if the formation of W UMa systems and their evolution 
(coalescence) into a single star is a continuous process, and 
somehow if there is an equilibrium in the Galaxy, which 
could be investigated of course, the 1.61 Gyr age difference 
between field CAB and field FCB implies an upper limit for 
the timescale of the contact stage. Indeed, it is not much 



different than the estimate of 0.1 < t conta ct < 1 Gyr by 
Guinan & Bradstreet (1988). 

Therefore, it is not possible to assume the MG group to 
be the initial stage of the FCB group. It is most likely that 
the MG contact systems are formed in the beginning of the 
main-sequence or during the pre-main-sequence contraction 
phase, either by a fission process (Roxburgh 1966) or most 
probably by fast spiraling in of two components in a common 
envelope. On the other hand, FCBs must be mostly formed 
from the detached progenitors. 

The evidence of total mass and orbital period decrease 
appears to be better displayed in Table 6 where the field 
contact binaries are divided into four sub groups according 
to orbital period criterion. It is clearly displayed there that 
the average kinematical age increases as the orbital periods 
of the sub groups decrease. It is also clear that older sub 
groups have smaller total mass. As if, long period contact 
binaries in Table 6 lose mass and angular momentum, so 
their orbital period decrease within the time scale indicated 
in the Table. This scenario, however, is definitely misleading 
and cannot be true. This is because the age difference be- 
tween the youngest and oldest sub groups in Table 6, is too 
long for a contact system to survive even if the empirical 
1.61 Gyr time scale for the contact stage were taken into 
consideration. 

According to the true scenario, the large dispersions 
(older ages) in Table 6 must have been produced not only 
during the contact stage, but also including the pre-contact 
stages before, which may be longer than the contact stage 
as indicated by Table 6. Thus, one should never expect, for 
a W UMa system in the youngest group in Table 6, to loose 
mass and angular momentum and join the shortest period 
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Table 5. Physical parameters of W UMas. 
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03 
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o 


52 
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1 
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o 
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Q 


70 
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1 1.25 


o 


44 


F5V 
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58.6 


1 


.176 
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(1 
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o 


22 


10 
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o 


62 
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75.3 


1 


.343 
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o 
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o 


13 


11 
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o 


78 
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.3050 


85 9 
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o 
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o 


03 
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7.68 


(1 


26 
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o 
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13 


EE Cct 


8.76 
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36 
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.3799 


77.5 
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.391 
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o 
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o 


20 


14 


UX Eri 
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Q 


67 
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.366 
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(oldest) W UMa systems after t = 8.89 - 3.21 = 5.68 Gyr 
later. It is not yet clear, but it seems that when a detached 
system becomes a contact binary, the angular momentum 
loss becomes accelerated with the gravitational radiation as 
described by Guinan & Bradstreet (1988). If the total mass 
of the system does not decrease to fit in the decreasing Roche 
lobes, then, overcontact / parameter starts to increase. The 
pre-contact stage, which may differ one system to system, 



being much longer than the contact stage, spoils the corre- 
lation between the / parameter and the space velocity as 
described in Fig. 5. 

Fig. 10 displays the total mass histograms of possible 
MG and FCB groups. Similar to the period histograms that 
it would be meaningless if one tries to deduce evidence of 
mass decrease during contact binary evolution from this fig- 
ure. Therefore, it should be considered as a bare comparison 
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Table 5 - continued 
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between MG and FCB. It is the same with Fig. 11, which 
compares the mass ratio distribution between MG and the 
FCB groups. The mass ratio evolution during the contact 
stage cannot be deduced from Fig. 11 since MG group can- 
not be the initial distribution. On the other hand, all four 
histograms (Figs. 8, 9, 10, and 11) display an empirical evi- 
dence of contact binary evolution to indicate coalescence of 
W UMa systems into single stars. Otherwise, angular mo- 



mentum loss does not fit in the theory of binary evolution. 
Then, it becomes difficult to explain why there are such dif- 
ferences between the MG and the FCB group W UMa stars. 

A similar case also applies to the Table 6 data. We be- 
lieve the age difference between MG contact binaries and 
the youngest FCB group in Table 6 is empirical evidence 
of coalescence into single stars. This is because the coales- 
cence into a single star involves orbital period decrease as 
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Figure 10. Comparison of the total mass (M^ + M c ) histograms 
of (a) MGs and (b) FCBs. 



Figure 11. Comparison of the mass ratio (q = M2/M1, where 
M2 < Mi, and Mi, M2 are the primary and secondary masses, 
respectively) histograms of (a) MGs and (b) FCBs. 



Table 6. Kincmatical ages for the period sub-groups in FCB. 



p 


< Mtot > 


< P > 


N 


Otot 


Age 


(days) 


(M ) 


(days) 




(km s _1 ) 


(Gyr) 


(0.5-0.9] 


2.220 


0.669 


17 


42.09 7.1 


3.21 


(0.4-0.5] 


1.900 


0.434 


34 


43.86 10.4 


3.51 


(0.3-0.4] 


1.618 


0.353 


35 


63.90 16.0 


7.14 


(0.2-0.3] 


1.171 


0.270 


11 


73.92 14.6 


8.89 



well as angular momentum loss, most probably accelerated 
by gravitational radiation as described by Guinan & Brad- 
street (1988) although the data in Table 6 is mostly affected 
by the long pre-contact stages. The derivation of the rates of 
orbital period decrease, mass and angular momentum loss 
from the present data of contact binaries requires a careful 
analysis and interpretation of Table 6 data. Due to limited 
space in this study derivation of those rates will be handled 
in a forthcoming study. 



5 CONCLUSION 

Kinematics of 129 W UMa binaries is studied. The sam- 
ple is found to be heterogeneous in the velocity space that 
kinematically younger and older contact binaries exist in 
the present sample. Various sub groups are formed accord- 
ing to criteria involving the spectral types, binary type (A 
or W), mass ratio, and over-contact parameter /, orbital 
period and total mass. However, those criteria failed to pro- 
duce kinematically homogeneous sub groups. Nevertheless, 
a correlation has been found to indicate shorter period and 



less massive systems show larger velocity dispersions than 
longer period more massive systems. The mean kinematical 
ages of all sub groups are estimated from the velocity disper- 
sions using the tables of Wielen (1982) and listed in Table 
3. 

The possible moving group members among the present 
W UMa stars are investigated according to the kinematical 
criteria originally defined by Eggcn (1958a, b, 1989, 1995) 
and 28 W UMa are found to be a possible member of five 
best documented MGs which are summarized in Table 4. 
The group of the field contact binaries (FCB) is formed by 
removing out the possible moving group members and the W 
UMa systems which are known to be the members of known 
open clusters from the main sample. The 5.47 Gyr mean 
kinematical age is estimated for the field contact systems 
(FCB). The age difference is found to be 1.61 Gyr between 
the 3.86 Gyr old CAB (Karatas, et al. 2004) and the FCB. 
If there is a number equilibrium between the CAB and FCB 
established by the evolution, the upper limit of the lifetime 
of the contact stage can be estimated as 1.61 Gyr which is 
in close agreement with the estimation made by Guinan & 
Bradstreet (1988) as 0.1 < t contac t < 1 Gyr. 

Since the contact stage appears to be very short in com- 
parison with the mean kinematical age of FCB, the group 
containing possible MG stars cannot be used as an initial sub 
sample of FCB. Thus, the histograms (Figs. 8, 9, and 10) and 
ages of sub groups of FCB according to orbital ranges (Ta- 
ble 6) cannot be used directly to deduce evidence of period 
decrease and angular momentum loss during contact stage 
evolution. Nevertheless, since, those histograms and data in 
Table 6, cannot be explained without angular momentum 
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loss and orbital period decrease in binary evolution, and 
since it is illogical to assume that angular momentum loss 
and orbital period decrease operate only for the pre-contact 
stages, the angular momentum loss and the orbital period 
decrease are evidently occuring as predicted by Guinan & 
Bradstreet (1988) even if one cannot deduce them directly 
from those histograms and Table 6. 

We have found 27 systems to form the MG group. Ex- 
istences of MG members is an exciting (perhaps the most 
exciting) finding of this study. This is because these W UMa 
systems do not have enough ages to be formed from detached 
binary progenitors by the mechanism involving angular mo- 
mentum loss and orbital period decrease. In fact the age 
difference between the youngest FCB sub group (3.21 Gyr, 
in Table 6) and MG group (0.5 Gyr, in Table 3) is more 
than lifetime (1.61 Gyr) of the contact stage. Therefore, it 
could be concluded that the MG contact systems are formed 
in the beginning of the main-sequence or during the pre- 
main-sequence contraction phase, either by a fission process 
(Roxburgh 1966) or by fast spiraling the two components in 
a common envelope. As a result, it can be said that detached 
binary progenitors, presumably short period (P < 5 days) 
RS CVn binaries, are not the only source from which W 
UMa systems could be formed by the mechanism involving 
angular momentum loss and orbital period decrease. It is 
indeed a challenge to disprove true membership of MG stars 
in our list. This may not be enough, one also needs to prove 
their ages to be bigger than our estimate, that is, their ages 
must be proven to be big enough (bigger than the pre-main 
sequence contraction plus the contact stage life time) then 
one can claim the detached binary progenitors are the only 
source to form W UMa systems. 
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